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1. Introduction 

Various low- dimensional electron systems have been the subject of considerable scientific 
interest for more than two decades. Examples of these structures can be formed in 
a two-dimensional electron gas by applying appropriate confinement, thus restricting 
the motion of electrons to a small area on the boundary between two semiconducting 
materials. These systems are in many ways similar to atoms, however, their properties 
can be controlled by adjusting their geometry, the external confinement, and applied 
magnetic field. Nanostructures are a source of discoveries of novel quantum phenomena 
which do not appear in atoms. They are important both in connection with potential 
device applications and can function as convenient samples to probe the properties of 
many-electron systems in reduced dimensions. 

In quantum rings, electrons are confined to move on a circular region in space. 
They have many interesting and actively explored properties: Under the infiuence 
of a magnetic field, an equilibrium current fiows along the ring. Furthermore, the 
energy spectrum is periodic in the magnetic fiux with a period of the fiux quantum 
$0 = h/e in case of an even number of noninteracting electrons, but $o/2 for an odd 
number. Recently, the existence of a fractional periodicity has been discovered both in 
experiments P as well as in theoretical and computational studies [21 El IH El 13 El ■ 
This is often referred to as the fractional Aharonov-Bohm effect. 

In this paper we study the accuracy of a variational quantum Monte Carlo 
(VMC) approach, using the Jastrow-Slater wave function, to quantum ring systems. 
This approach has been found to be extremely accurate in the case of quantum 
dots [9 . The main advantages of quantum Monte Carlo (QMC) methods are their 
applicability to large electron systems and their accuracy in capturing electron-electron 
correlation effects. To date, very few QMC studies of the electronic properties of 
nanoscale semiconductor quantum rings have been done. Emperador et al investigated 
the fractional Aharonov-Bohm effect in few-electron quantum rings using multi- 
configurational diffusion Monte Carlo [21 and Borrmann and Harting used path integral 
Monte Carlo to study the transition between spin ordered and disordered Wigner crystal 
rings |10j. We are not aware of VMC calculations for the quantum rings, and due to this, 
even the accuracy of the Jastrow-Slater wave function is an interesting open question. 
On the other hand, the good performance of VMC on two-dimensional quantum dots 
makes it an appealing technique to study also systems that are closer to being one 
dimensional fH^l. As an ultimate test for VMC in this direction, we consider here 
the limit of a purely one-dimensional system. 

We calculate using VMC the total energy and pair distribution function of the 
quantum ring, in zero magnetic field, containing up to six interacting electrons. The 
interaction is taken to be Coulombic. To test the accuracy of our approach we perform 
exact diagonalization (ED) calculations and compare the results. We will show that 
VMC gives a very accurate estimate of the ground state energy for all the cases studied 
here. In addition, the charge part of electron correlations is recovered in very satisfactory 
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fashion. The only weak point of VMC is that the antiferromagnetic nature of the system 
seen in the exact treatment is not present in the VMC resuhs. This shows that even 
though the simple VMC wave function is adequate for the charge part, the spin structure 
of the Jastrow-Slater wave function is not correct. On the other hand, the spin part has 
only a minor role in, e.g., the total energy. 

Already at a moderate strength of interaction, the few-electron quantum ring shows 
a separation of spin and charge degrees of freedom. The spin Hamiltonian is well 
described by the antiferromagnetic Heisenberg model with exchange coupling J jl3j . 
which is very small compared to the energy of orbital motion. We compare our results 
to the relations for J derived by Fogler and Pivovarov |E] and show that they are in 
good agreement. We discuss the localization of electrons in relative coordinates and find 
a smooth transition starting already at 1 < < 2. 

This paper is divided as follows. In section II we discuss our model of the quantum 
ring system and go in considerable detail through the two methods used to calculate its 
ground state properties, ED and VMC. Then the results are presented in Section III. 
Our conclusions follow in Section IV. 

2. Model and methods 

2.1. Model 

We model the semiconductor quantum rings in the effective mass approximation by a 
many-body Hamiltonian 



where 6i is the angular coordinate of the ith electron and V{6ij) is the interaction 
potential between each pair of electrons. We assume that the ring is extremely narrow so 
that the electrons can only move at a certain radius, leading to a purely one- dimensional 
ring. The Coulomb interaction for a ring system is written as 



where R is the radius of the ring, Vc controls the strength of the interaction, and /i is a 
small parameter that eliminates the singularity at fty = 0. As the interaction between 
particles at the same position on the ring is Vc/fiR, one can think the regularization of 
the potential to give the ring a width of fiR. 

The single-particle eigenstates of the noninteracting case are plane waves 




(1) 



Vc 



(2) 



R^/Asin' (%/2)+/i2 ' 



(3) 



with energies 



£i = 



2m*R^ 




(4) 
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where / = 0, ±1, ±2, ... is the angular momentum quantum number. We scale the units 
so that the energy is measured in units of Eq = /i^/(2m*7r^i?^) and length in units of 
R. The only parameter which is tuned is the interaction strength, and for the Coulomb 
interaction this can be translated into a change of the ring radius by the formula 

where is the effective Bohr radius. The results presented in this paper are for 
interaction strength in the range of 1 to 10, which corresponds to a radius of 21 nm to 
209 nm assuming GaAs material parameters. For comparison, the quantum ring system 
studied by Keyser et al ^ had a radius of around 150 nm, which corresponds to an 
interaction strength of Vc = 7.2 in our units. The width of their ring can be estimated 
to be around 20 nm, which is in reasonable agreement with the value of /i = 0.1 used 
in this work. 

In the following, we will occasionally use the one-dimensional r^-parameter, which 
is the average distance between electrons. In our model, it is directly proportional to 
the radius of the ring and thus also to Vc 



a*B N 



(6) 



2.2. Exact diagonalization 

The exact diagonalization method (ED), originally called configuration interaction, is a 
systematic scheme to expand the many-particle wave function using the noninteracting 
states as a basis. This method traces back to the early days of quantum mechanics, 
to the work of Hylleraas ^Hl- The calculation of matrix elements for the many-body 
case was originally derived by Slater and Condon [T71 UHl HH] , and developed further by 
Lowdin [20] • The use of the term ED for this can be seen as an attempt to replace the 
term configuration interaction with a normal mathematical term [25 . This term might 
in some cases be misleading, as truly exact results are obtained only in the limit of an 
infinite basis. 

As a simple example to start from, consider a single-particle Hamiltonian split into 
two parts Ti = Hq + Hi, where the Schrodinger equation of the first part is solvable: 

7^o0^(r) = £.0^(r) , (7) 

and the wave functions 0i(r) form a orthonormal basis. The solution of the full 
Schrodinger equation can be expanded in this basis as ^^(r) = ^^aj0j(r). Inserting 
this into the Schrodinger equation 

n^{r) = Eij{r) , (8) 

results in 

(Ho + Hi) J2 = eY. "^'^^(^) • (9) 
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Using (I7|), multiplying with 0* and integrating gives 

^OiEidij + ^ai / (j)*{r)Hi(f)i{r)dr = E^aiS, 



(10) 



i i i 



for every j. This can be written as a matrix equation 
{Hq + Hi)cx = Ea , 



(11) 



where Hq is a diagonal matrix whose ith element is Ei, and the jith element of Hj is 



equation can be mapped to a matrix form. In principle, the basis {(pi} is infinite, but 
the calculations are done in a finite basis. The main computational task is to calculate 
the matrix elements of Hj and to diagonalize the matrices. The convergence of the 
expansion depends on the actual values of the matrix elements, and in the case where 
Til is only a small perturbation to TYq it is fastest. 

In a similar fashion, one can split the many-particle Hamiltonian into two parts, 
typically a noninteracting and an interacting one. The solution to the noninteracting 
many-Fermion problem is a Slater determinant formed from the N lowest energy 
eigenstates of the single-particle Hamiltonian. The basis typically chosen for the 
interacting problem is a Slater determinant basis constructed from M single-particle 
states. For an A^-particle problem, ('^) different many-particle configurations can be 
constructed by occupying any of the M available states. Each of these determinants 
is labeled by the N indices of the occupied single-particle states. Due to interactions, 
other configurations than the one of the noninteracting ground state have a finite weight 
in the expansion of the many-particle wave function. 

The Schrodinger equation maps to matrix form in a similar way as in the preliminary 
example above. The calculation of the matrix elements of the Hamiltonian are simpler 
using second quantization and the occupation number representation. We will only state 
the results here. 

The noninteracting Hamiltonian matrix Hq in the determinant basis is diagonal 
with the zth element equal to Ei = J2f=i ^ij 5 where ij labels the occupied single-particle 
state. The elements of the interaction Hamiltonian matrix can be found by a straight- 
forward calculation using the anti-commutation rules of the creation and annihilation 
operators. The only non-zero matrix elements are between configurations that differ 
at most by two single-particle occupations. If both of the two configurations have two 
single-particle states occupied that are unoccupied in the other one, the interaction 
matrix element is 



where p and q are single-particle states occupied in the left configuration and unoccupied 
in the right one, and similarly for s and r with right and left interchanged. The sign of 
the matrix element depends on the total number of occupations between p and r and 




contains the values ctj. In this way, the Schrodinger 




(12) 
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between q and s due to anti-commutation of the creation and annihilation operators. If 
this number is odd, the sign is minus, otherwise plus. For the Coulomb interaction 

Vijkl = Ss,.sjs,,si I 0*(ri)0*(r2)— 0fc(ri)0;(r2)dridr2 , (13) 

J ^12 

and the delta functions result from the spin part summation. In the case of a one- 
dimensional quantum ring this matrix element has a closed form 

Vjki = (/iV2 + 1) , (14) 

where Q{x) is the Legendre function which can be written in terms of Gauss's 
hypergeometric function. 

If the two configurations have a difference of only one occupation, the final result 
for the interaction matrix elements is 

i ^^{Vqiri Vqiir) , (15) 



where q is occupied in the left configuration and unoccupied in the right one, and for 
r the other way around. The sum over i is over orbitals that are occupied in both 
configurations. Again the sign of the matrix element depends on number of occupations 
between the differing orbitals. It is easy to see that the matrix elements in (fT3j) vanish 
for a rotationally symmetric system, as a change in only one occupation cannot conserve 
the total angular momentum. 

The diagonal elements of the interaction Hamiltonian matrix are simply 

N 
i<j 

To summarize, the noninteracting part results in a Hamiltonian matrix which is 
diagonal, and the interaction Hamiltonian couples configurations that differ by at most 
two occupations. Now that we know the Hamiltonian matrix elements, we can sketch 
the basic ED procedure as follows: First solve for M eigenstates of the single-particle 
Hamiltonian Tio- After those have been found, calculate the two-particle matrix elements 
Vijkl of (|T!?|) . Construct the A^-particle configurations from M single-particle states. 
Here configurations with wrong symmetry can be rejected. For example, one can pick 
only configurations with certain z-component of the total spin or some other good 
quantum number, like the angular momentum in our case. To construct the Hamiltonian 
matrix, calculate the interaction matrix elements of the Hamiltonian Tij between the 
configurations. Construct the diagonal noninteracting Hamiltonian matrix elements 
from the single-particle eigenvalues. Finally, diagonalize the Hamiltonian matrix. The 
main problem of the ED method is the exponential computational scaling of the method 
as a function of the number of particles. An other problem is the convergence rate as a 
function of basis size. 

To compare the ED and VMC methods for the ring, we study the total energy and 
the pair distribution function of the electrons in the ring. The pair distribution function 
of the system is defined by 

n^(r)n.,(r')w(r,r') = (^| ^i{r)^i,{r')^^,{r')^^{r) |^) , (17) 
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where no-(r) is the spin a particle density and \&J.(r) and \I/cr(r) are the field operators. 
It describes the probability that, given an electron of spin a at r, an electron of spin a' 
is found at r'. It therefore gives an idea of the strength of correlation between electrons. 
The pair distribution function can be derived in a similar way as the Coulomb matrix 
elements using the ant i- commutation rules of the creation and annihilation operators. 



2.3. Variational Monte Carlo method 

The QMC methods are among the most accurate ones for tackling a problem of 
interacting quantum particles [22j. Often the simplest QMC method, namely, the 
variational QMC (VMC), is able to reveal the most important correlation effects. In 
many quantum systems, further accuracy in, e.g., the energy is needed, and in these 
cases methods such as the diffusion QMC allow one to obtain more accurate estimates 
for various observables. 

The VMC strategy of solving the ground state of an interacting system of N 
particles (defined by the Hamiltonian operator H and the particle statistics) starts 
by constructing the variational many-body wave function \&. Then, all the observables 
of the system can be computed from it. For example, the energy E, which is higher 
than the ground state energy Eq, can be obtained from the high- dimensional integrals 

/|M/(R)|MR ' ^^^^ 
where R is a vector containing all the coordinates of the particles. If the particles 
have d degrees of freedom, the integrals here are N x d dimensional. As one would like 
to have reasonable large, the integrals are typically so high-dimensional that these 
are most efficiently calculated using a Monte Carlo strategy. This is based on the limit 

Nr^oo A^R ^ V Jv 

where the points Rj are uniformly distributed throughout the volume V. The error in 
the evaluation of the integral decays as oc o-f/y/Nji, where cx/ is the standard deviation 
of the function /. A typical trick to reduce the error is to make ctj smaller by dividing 
the function / by a similar function g (that we assume to integrate to one). Then, one 
can rewrite the sum as 

^|^/(R,)-^1^9(R.)j^-^l^^, (20) 

where now the points Rf are distributed as g. The error in this approach is now smaller, 
as the function g was chosen such that the standard deviation of f/g is smaller than 
that of /. 

Now in VMC, one typically generates the random set of Nji coordinates {Rj}^^ so 
that they are distributed according to |\E'p, and one obtains for the energy 

1 ^ nim _ 1 ^ 
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(22) 



^(R.) 



One should note that El is a function of the coordinates, and it defines the energy at this 
point in space. If the wave function \l/ solves the Schrodinger equation, El equals the 
true energy of the system. From the Monte Carlo point of view, the error in determining 
the energy is related to the standard deviation of the local energy, and for this reason 
the statistical error in the VMC energy decreases as the wave function becomes more 
accurate. 

The sampling of the random points Rj can be done by, e.g., using the basic 
Metropolis algorithm. In that, the ratio of the sampling function at the ith and the 
previous step: |\E^(Rj)p/|^(Rj_i)p, needs to be calculated and the trial step is accepted 
if the ratio is larger than a uniformly distributed random number between zero and one. 

A typical, and also the most commonly used VMC trial many-body wave function, 
is the one of a Slater- Jastrow type: 



where the two first factors are Slater determinants for the two spin types, and 



where is a Jastrow two-body correlation factor. The determinants contain 

single-particle orbitals, where the ones in different spin determinants can be the same. 
These determinants solve in some cases a noninteracting or a mean-field Schrodinger 
equation, but in general, one can choose these orbitals rather freely. This form of a 
wave function has proved to be very accurate in many cases 12^. In quantum dots, the 
Jastrow- Slater wave function captures most correlation effects very accurately [Hj. The 
only possible exceptions are found in the fractional quantum Hall effect regime, where 
only some of the states can be well approximated by the Jastrow correlation. Interesting 
examples of successful cases are the Laughlin states P^j . 

One can generalize the Jastrow-Slater wave function to contain more than one 
determinant per spin type. This leads to an increased computational complexity, 
and for this reason it is often avoided. Another generalization can be to replace the 
single-particle coordinates in the determinants by a set of collective coordinates, e.g., 
Ti ^ Ti — F{{Tj}^^j). This also leads to complications, this time in the calculation of 
the ratio of the wave functions in the Metropolis algorithm. 

The determinants in our VMC trial wave function are constructed from the 
noninteracting single-particle states, and for the two-body Jastrow factor we use here a 
simple form of 



n 



(23) 




(24) 



i<j 



10 




(25) 



i=0 
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where a's are variational parameters, different for electron pairs of same and opposite 
spin type. 

An important ingredient in VMC is the optimization of the variational parameters 
of the wave function |j9j . For an efficient optimization of the variational parameters, we 
need the derivative of the energy with respect to these parameters. Our wave function 
is complex-valued, so we need to generalize the result for the real case [21] to read: 

dE d J **7^^de 



dak 



dak /|^|M0 
1 



/ |^|2de 
1 



(/|*|Me)2 

2 



(/|^|2de)2 



^*7^^de 



de 

V]/*V]/' -)- ( \[/*v[/ 



de 



— ^de 



1^ 



1^1 



de 



2SR 



E, 



(26) 



where = d'^/dak and de = d^^i, . . . , d^Ar. In the second line of the equation we 
have used the fact that the Hamiltonian Ti is Hermitian and real. The angle brackets 
denote an average over the probability distribution l^&p. This notation is also used 
in the following. In our wave function, the parameters ai are inside an exponential, 
real-valued Jastrow factor: 

where e = (6'i, . . . , 9^), a = (ai, . . . , q;m )• As a result we have: 

^ dak 
Using PH|) we can rewrite as 
dE 



(27) 



f28) 



dai 



23? 



Q Er 



Q )(E 



(29) 



3. Results 



3.1. Convergence and configurations 

Before presenting the actual results we first present an analysis of the accuracy of the 
ED method used. The only approximation done in our ED calculations is that we have 
a finite number of determinants in our expansion of the many-electron wave function. 
The actual number used is finally limited by the available computer resources. To give 
an idea of the accuracy of ED, we show in figure Q the convergence of the = 6, = 



Interacting electrons on a quantum ring: exact and variational approach 



10 



ground state energy with increasing basis size M. The interaction strength is = 10 
making this the worst-case scenario for the convergence. We have increased M in steps 
of four, because the energy is lowered most when a new pair of ±/ angular momentum 
values are included for both spin types. In this case, the maximum number of basis 
functions we have used is M = 46. For these values of M and N, {^^) = 9 366 819. 
However, many configurations can be rejected from symmetry requirements, as we can 
force Sz and L to have certain fixed values. This use of symmetry leads to a final size 
of the Hamiltonian matrix of around 79 000 x 79 000. The difference in energy between 
M = 42 and M = 46 is 4.2 ■ 10"^ (in units of Eq) which is about 0.00034% of the total 
energy. We can therefore be confident that our ED calculations have sufficient accuracy 
to be used to test the quality of the VMC approach. 

When the electrons interact weakly, the many body state has one major 
configuration, that of the equivalent noninteracting system. In this limit, the single- 
determinant Hartree-Fock (HF) method gives reasonable results. The electrons arrange 
so as to minimize the total kinetic energy and the single-particle states are occupied as 
compactly as possible around the / = state. The resulting configurations are shown in 
the upper part of figure El Each box in this figure represents one single-particle state, 
having angular momentum from / = at the bottom, up through / = ±1, ±2, . . ., with 
increasing energy. For four electrons, the second shell formed by the angular momentum 
/ = ±1 states is half-filled, and the exchange energy favours spin-polarization of the two 
electrons in that shell. When the interaction strength grows, the relative importance of 
the kinetic energy decreases with respect to the interaction energy. The most important 
configuration of the many-body wavefunction of strongly interacting electrons has closed 
shells, and thus zero angular momentum, for each spin type. The configurations are then 
either a half-filled shell of one spin type or a closed shell. The resulting configurations 
are shown in the lower half of figure El For = 4 and 6 electrons, the state is already 
a closed-shell at weak interaction and it remains the major configuration in the limit 
of strong interaction. The odd cases have open-shells at weak interaction and make 
a transition to closed-shell configuration at some point. For a three-electron ring, ED 
calculations show this transition to occur at an interaction strength of about Vc ~ 0.41, 
corresponding to ^ 0.85 a^. For five electrons the transition takes place at Vc ~ 3.3 
{rg ~ 4.1 a^). The determinant part of the VMC wavefunction is chosen to be the 
relevant configuration shown in figure |21 These are the main configurations of the 
many-body state found by ED and have, at all except the strongest interaction, by far 
the largest weight in the determinant expansion. 

3.2. VMC accuracy 

The total energy of six electrons on the ring is shown in figure El We have done the 
calculations for the L = state, and a total spin S = (groundstate) and S = 2, 
and show the results obtained from ED, VMC and HF as a function of the strength 
of interaction (proportional to the ring radius). Over the whole range, ED and VMC 
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Figure 1. Convergence of the ground state energy as a function of the number of 
single-particle states M, for the case of iV = 6 electrons and a strong interaction 
Vc — 10. The dashed line is to guide the eye. 
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Figure 2. Single-particle states occupied in the VMC wave function. These are also 
the most important configurations in the ED expansion. 



results agree very well, although the agreement improves slightly with increasing Vc- 
The Hartree-Fock results are reasonably accurate only at the lowest interaction strength 
and already at Vc = 2 (corresponding to R = Aa*^ and = 2.1 a^) shows considerable 
deviation from the ED results. 

To compare in more detail the results of ED and VMC, a convenient measure of 
the quality of the trial wave function is the percentage of correlation energy captured 
by VMC, defined as 

Ec = ^^^i^ , (30) 

i^HF — -C^ED 

where E^y is the Hartree-Fock energy. In figure 0] we show the results for the L = 
state of six electrons. At all interaction strengths, VMC captures most of the correlation 
energy. It is more accurate for the S = 2 state than the S = state, although the 
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Figure 3. Total energy as a function of interaction strength for an = 6 electron 
ring of spin S — Q and S — 1 and angular momentum L = 0. Results from ED, VMC 
and HF are shown. 
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Figure 4. Percentage of correlation energy captured by VMC [see (|30|l for definition] 
for = 6 electrons, with L = 0. The error in the the first two points is given by the 
marker size, and is smaller than that for the remaining points. 



difference decreases with increasing strength of interaction. One would expect VMC 
to become less accurate for strong interactions, because the construction of the VMC 
wave function starts from a single configuration which is later multiplied by a Jastrow 
factor, whereas the many-body wavefunction contains multiple determinants and their 
weights get more evenly distributed with increasing Vc- This is not the case, and 
at the highest Vc the deviation in the VMC groundstate energy lacks such a single- 
configuration limitation. For other particle numbers the results are shown in table 
for low and high interaction strength, Vc = 1 and Vc = 10. The VMC results are 
very accurate for the particle numbers we have considered. In the worst case, VMC 
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Figure 5. Pair distribution function from ED and VMC for = 6 electrons and 
L = with 5 = in gray and S = 2 in black. On the left Vc ~ I and on the right 
Vc — 10. The dots show the VMC data while solid lines are the ED results. For strong 
interaction, the curves are nearly identical. 



Table 1. Percentages of correlation energy captured by VMC. The error in the last 
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captures 96% of the correlation energy and the total energy is roughly 0.7% higher 
than that found by ED. The accuracy improves with increasing interaction strength 
and is better for lower particle numbers. This trend is broken by cases where the major 
configuration has an open shell, which is always the case for odd particle numbers at 
weak interaction. This is demonstrated by five electrons, where VMC captures slightly 
less of the correlation energy than for six electrons. 

In figure El the total pair distribution function of an = 6 electron ring is 
shown at weak and strong interaction. As in the case of the groundstate energy, the 
agreement between VMC and ED is good. The left panel shows clearly how the total 
pair distribution functions of different spin states become identical when the electrons 
interact strongly. This is one indication of spin-charge separation, when the spin and 
charge of a one- dimensional system decouple. 

The spin channels in the pair distribution function are shown in figure IHl For the 
S = 2 state, VMC results are almost identical to those from ED. In the case of S" = 0, 
however, the internal structure is inaccurate at weak interaction and completely wrong 
at strong interaction. This is further shown in figure [7| where we have plotted 

Ag = -^-^-^, (31) 

which for an antiferromagnetic alignment of spins should give equidistant peaks (up 
spins) and troughs (down spins). This is indeed the case, and we furthermore see a 
decay in spin correlations with increasing interelectron distance. The VMC results show 
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Figure 6. Spin cliannels in the pair distribution function from ED and VMC for 
= 6 electrons and i = 0. The upper panel shows the results for the 5 = state 
and interaction strengths Vc = 1 (left) and Vc = 10 (right). The lower panel shows 
results for the 5 = 2 state Vc = 1 (left) and Vc = 10 (right). The markers show the 
VMC data while solid lines are the ED results. 



no antiferromagnetic structure. We believe the reason for VMC not being able to handle 
antiferromagnetic spin coupling is due to three-particle correlation which is missing in 
the VMC wavefunction. Despite VMC's shortcomings in accurately describing the spin 
correlations for an antiferromagnetic spin structure, the total pair distribution function 
and the groundstate energy are correctly described. 

3.3. Spin and charge in the Wigner crystal regime 

Like the pair distribution function, the total energy of the S = and 5* = 2 states 
becomes identical in the limit of strong interaction. In figure |H1 we show the energy 
difference of these states for six electrons, as a function of Vc found by ED, VMC and 
HF. The inset shows the energy difference found by ED on a logarithmic scale. As 
expected, HF fails seriously already at low values of Vc and predicts a groundstate of 
nonzero spin-polarization at Vc just below 1 (r^ ^ 1). VMC follows the ED results 
more closely, but at Vc- = 4 predicts a change in spin. ED results indicate, however, 
that the groundstate spin remains S = over the whole range of Vc although the 
energy difference approaches zero. Interestingly, the difference between the ED and 
VMC results remains approximately constant for Vc- = 4 . . . 10. 
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Figure 7. Pair distribution function difference, Ag — g|f/(A^| — 1) — gn/Ni, 
for six electrons in the state L — 0, S ~ 0. The ED results clearly indicate an 
antiferromagnetic structure. 



Figure |H1 shows that in the strong interaction regime, the energy scale of spin 
dynamics becomes very small when compared to the energy scale of the orbital motion. 
This is one indication of a strong spin-charge separation in the system. It is further 
manifested in the pair distribution function, as shown in figure which is the same 
for both spin states, S = and S' = 2, at stronger interaction. The spin part can be 
described by the antiferromagnetic Heisenberg model ^3] with spin interaction 

= J^Si-Sj (32) 

with a nearest-neighbour coupling constant J. In a recent paper, Fogler and Pivovarov 
derived a relation between the coupling constant of thin quantum rings and the pair 
distribution function of a spin-polarized ID system go{x) 

J = g''(o) (33) 

where R* is the thickness of the ring, in our case /ii?. For ultrathin wires exhibiting 
strong spin-charge separation, this relation holds even at moderate electron densities, 
~ 1. They furthermore derive an expression for the pair distribution function, which 
leads to 

K TT e^ 

= /o ^SM W exp(-r7V2^) , (34) 

where k and rj are constants. The coupling constant J is proportional to the difference 
in energy of the S = and S = 2 states, AE. In our model, = 2ttVc/N so we can 
write 

AE oc exp(-r7v^47r/A^yVb) • (35) 
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Figure 8. Difference in energy of the S* = and S* = 2 states of six electrons with 
angular momentum L = 0, as found by ED, VMC and HF. In the inset the logarithm 
of the energy difference is plotted against VVc- 
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Figure 9. Pair distribution function of six electrons in the L — 0, S = state for 

several values of the strength of interaction Ve- 
in the inset of figure |H1 we have plotted the logarithm of AE as a function of a/Vc'- Our 
results agree very well with the results of Fogler et al., as the curve becomes linear even 
at low values of Vc, Vq ~ 1-7 (r^ ~ 1.8). From our data we can estimate the slope, 
and find that ?7 ~ 2.5 which is reasonably close to the value 2.7979 found by Fogler and 
Pivarov. However, it must be noted that our results are for a fixed value of [i = 0.1, 
so that the thickness of the ring scales linearly with increasing interaction strength as 
R*/a% = O.IR = 0.2Vc. 

The pair distribution function of six electrons of zero total angular momentum 
and spin is shown in figure El for several values of Vc- At Vc- = 1 the distribution is 
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Figure 10. Logarithm of the contact probabihty gniO) of opposite spin electrons as a 
function of the strength of interaction Vc on a square root scale. The parabolic curve 
is a guide to the eye only and indicates an exponential fallofF with Vc- The straight 
line is a linear fit and indicates a exp(— ^-^ Vc) behaviour. 

rather flat, but already at Vc = 2 there is a clear peak structure indicating partial 
localization in relative coordinates. As the system is finite, a phase transition to a 
Wigner crystal does not take place. However, there is a smooth transition to a regime 
where a Wigner-crystal-like behaviour is found. The system is translationally invariant 
and the electron density remains uniform, the localization occurs in the relative frame. 
In two-dimensional quantum dots, the electrons are found to be localized at |25j = 4, 
and the transition is preceded by a change in the groundstate spin j^HI- A QMC study 
of the three-dimensional electron gas showed an abrupt change in the evolution of the 
correlation energy with density [2ZI- We do not see such a change in our system. One 
should also note that the three-dimensional QMC work has different wave functions for 
the liquid and sohd phase, whereas we use the same wave function in both cases. 

Further insight is obtained from the value of the pair distribution function at ^ = 0, 
which is the contact probability. At low values of the strength of interaction Vc it has 
a nonzero value, as electrons of opposite spin are not restricted by the Pauli exclusion 
principle. As Vc grows, the strong Coulomb repulsion reduces g^iiO) to the point where 
the electrons become localized between their neighbors on either side. A closer look at 
the contact probability reveals an indication that the quantum ring system undergoes 
a smooth transition at densities between 1 < < 2. In figure El the logarithm of the 
contact probability is shown as a function of the interaction strength on a square root 
scale. For Vc > 2, the linear fit indicates that 5'|j(0) falls off as exp(— Ca/T^), where ( 
is a constant. At low interaction strength, the behaviour is closer to exp{—('Vc), the 
parabolic curve shown. A related behaviour has been found for a different quantity in 
quantum dots 
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4. Conclusions 

We have studied interacting electrons in a one-dimensional quantum ring. We compared 
the results obtained from variational Monte Carlo with those of exact diagonalization. 
The variational wave function was taken to be of a simple Slater- Jastrow form: a single 
determinant, the ground state of the noninteracting system, multiplied by a two-body 
correlation factor. The results obtained indicate that this simple wave function is able to 
capture the important correlation effects, and in most cases even with a high accuracy. 

The major source of error in VMC is the spin energy scale and spin correlation. 
The spin energy scale, however, becomes very small already at rather weak interaction 
and the error in the groundstate energy is correspondingly small and decreases with 
increasing Vc- For an antiferromagnetic spin correlation, VMC does not capture the 
internal spin structure, but the total pair distribution function is accurate. From this, 
one can conclude that the Jastrow factor is very efficient in capturing the total energy 
and essential correlation effects induced by the interaction between electrons. One might 
question the importance of the determinant part of the VMC wave function. We believe 
that having a good choice for this is crucial for the success of the VMC simulation. 
Even if the weight of this single configuration is not large in the exact expansion when 
the interaction is strong, it specifies the symmetry of the VMC many-body state. One 
should note that this is also the case for quantum dots, where the noninteracting single- 
particle states can be used in the determinant part of the VMC wave function, and they 
are even the optimal ones in the case of a closed-shell configuration [H] . The reason, both 
in quantum rings and dots, is the high symmetry of the problem: In a parabolic dot, 
the centre-of-mass and relative motion decouple, also in the many-body wave function. 
Thus if one starts from a noninteracting state that is a ground state of the centre-of- 
mass motion, the interactions only change the relative motion part of the many-body 
wave function. This means that if a Slater- Jastrow wave function is used, the Slater 
determinant part is not changed by the interaction between electrons. The case of a 
quantum ring is similar in the sense that the centre-of-mass and relative motion also 
decouple, but now because of the translational invariance of the system. 

For a moderate strength of interaction, the spin correlation is well described by 
an antiferromagnetic Heisenberg model with coupling constant J. We showed that our 
results agree well with the expression for J derived by Fogler and Pivovarov ^3] at 
Ts ^ 2. Their result is applicable already at ~ 1 if a strong spin-charge separation is 
present in the system. This indicates that for an ultrathin quantum ring this separation 
does take place at a quite high electron density. The contact probability also shows that 
there is a change in the properties of the system starting at 1 < < 2. 

In conclusion, a simple variational wave function, combined with Monte Carlo 
techniques to calculate the optimal wave function and observables of the system, results 
in an efficient computational strategy for interacting electrons in purely one- dimensional 
quantum rings. Since a similar conclusion can be drawn for two-dimensional quantum 
dots ^9j, we expect this to hold also for quantum rings of a finite width. 
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